Magnetization distribution in the transverse Ising chain with energy flux 



V. Eislei-i, Z. Raczi^2 ^nd F. van Wijland^ 
^Institute for Theoretical Physics, Edtvos University, 1117 Budapest, Pdzmdny setdny 1/a, Hungary 
^ Laboratoire de Physique Theorique, Bat. 210, Universite de Paris-Sud, 91405 Orsay Cedex, France 

(Dated: January 23, 2003) 

The zero-temperature transverse Ising chain carrying an energy flux Je is studied with the aim 
of determining the nonequihbrium distribution functions, P{Mz) and P{Mx), of its transverse and 
longitudinal magnetizations, respectively. An exact calculation reveals that P{Mz) is a Gaussian 
both at Je ~ and Je 0, and the width of the distribution decreases with increasing energy 
flux. The distribution of the order-parameter fluctuations, P{Mx), is evaluated numerically for 
spin-chains of up to 20 spins. For the equilibrium case {Je = 0), we find the expected Gaussian 
fluctuations away from the critical point while the critical order-parameter fluctuations are shown to 
be non-gaussian with a scaling function $(a;) = ^{M^/ {Mx)) ~ {Mx)P{Mx) strongly dependent on 
the boundary conditions. When Je 0, the system displays long-range, oscillating correlations but 
P{Mx) is a Gaussian nevertheless, and the width of the Gaussian decreases with increasing je- In 
particular, we find that, at critical transverse field, the width has a j'^^^^ asymptotic in the Je 
limit. 
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I. INTRODUCTION 

Nonequilibrium steady states (NESS) have been much 
studied but a description of some generality has not 
emerged so far. Among the many approaches tried, 
there are two which continue to receive particular at- 
tention. One of them is an attempt to understand the 
general features of NESS through studies of nonequilib- 
rium phase transitions 0, Qi 01 ■ The basic assumption 
here is that the universality displayed by equilibrium 
phase transitions carries over to critical phenomena in 
NESS, as well. Thus, by investigating the similarities 
and differences from equilibrium, one may gain an under- 
standing of the role of various components of the com- 
peting dynamics generating the steady state. For exam- 
ple, one may find by observing the universality classes of 
various nonequilibrium phase transitions that dynamical 
anisotropics often yield dipolc-likc effective interactions 
0, [5I, [6| or that competing non-local dynamics (anoma- 
lous diffusion) generates long-range, power-law interac- 
tions 0. 

The second approach is more ambitious but less well 
founded. It concerns the building (approximating, guess- 
ing) of the entropy or the distribution function for the 
nonequilibrium microstates 0, 0, . The idea here is 
that NESS is often the result of complex nonlinear dy- 
namics which yield distributions that cannot be described 
by the Boltzmann-Gibbs statistics. The reason e.g. may 
be that the dynamics generates long-range correlations 
which results in the loss of additivity of such quantities 
as the entropy. Furthermore, the dynamics may con- 
centrate the distribution onto a fractal subspace in the 
long-time limit, a case that again may be problematic in 
equilibrium statistics. In order to treat such cases, the 
so called non-extensive statistical mechanics was intro- 
duced 1^]. It is an approach that takes its name from 
the non-extensive character of the postulated entropy. It 



has been much developed [31 as well as criticized (see e.g. 
[iTl IT^ for the latest examples) during the last decade. 
Not surprisingly, the approach had its success in connec- 
tion with systems which have long-range interactions or 
display (multi)fractal behavior P|. 

Our work presented below is perhaps best viewed as an 
attempt to build a bridge between the above approaches. 
On the one hand, we study the effects of a nonequilib- 
rium constraint on a well investigated phase transition. 
Namely, we take the transverse Ising chain which has 
an order-disorder transition as the transverse field {h) is 
varied, and drive it by a field to produce an energy flux, 
Je through the system. The resulting steady states have 
been investigated ^T§\ and it has been found that, in addi- 
tion to the equilibrium phases, a flux-carrying nonequi- 
librium phase appears which is distinct by its correla- 
tions decaying with distance as a power law. On the 
other hand, we shall be concerned with the distribution 
function in the various phases of the above system. More 
precisely, we shall determine the steady-state distribution 
functions, P{Mz) and P{Mx) of the Mz (non-ordering 
field) and (ordering field) components of the macro- 
scopic magnetization in all three phases of the system 
and at and near its critical point. 

The results are surprisingly simple. The distribution 
functions are Gaussian in the equilibrium phases away 
from the critical point. This is expected since we have 
macroscopic quantity and the correlations decay expo- 
nentially. The distribution of the non-ordering field re- 
mains a Gaussian at the critical point of the equilibrium 
system, as well. The reason for this is that although 
the appropriate correlations decay with distance n as a 
power law but the exponent in the power is large (1/n^), 
so that the fluctuations (M^) — (M^)^ do not diverge at 
h = he- The distribution function of the ordering field 
becomes nontrivial at he and our numerical calculations 
demonstrate that P{Mx) depends strongly on the bound- 



2 



ary conditions taken to be periodic, anti-periodic, and 
free. The unexpected simplicity is in the current-carrying 
phase where the energy flux generates long-range correla- 
tions decaying as a power law {l/^/n) but, nevertheless, 
the distributions are gaussian. The mathematical reason 
for this lies in the oscillating character of the correlations 
which prevents the divergence of the spatial sum of the 
correlations which in turn are proportional to the fluctu- 
ations. Physically, the oscillations in the correlations can 
be traced to the form of energy flux [see eq.© below] 
which suggest that the consecutive x and y components 
of the spins are more and more rigidly interconnected 
as jE is increased and thus fluctuations decrease with in- 
creasing jE- This picture will be seen to be valid near the 
nonequilibrium phase boundaries where the fluctuations 
as a function of je can be explicitely calculated. 

The above results are presented in the following or- 
der. Section Hll contains a review of the transverse Ising 
model with energy flux, including the setup of the formal- 
ism convenient for calculating the distribution functions. 
Next (Sec lIII|l . the distribution P{Mz) is calculated ex- 
actly. Numerical work on P{Mx) and preliminary ana- 
lytic work on correlations are presented in Sec II VI fol- 
lowed by concluding remarks in SeclVI 



II. TRANSVERSE ISING MODEL WITH 
ENERGY FLUX 

The transverse Ising chain is one of the simplest system 
displaying a critical order-disorder transition )15j. It is 
defined by the Hamiltonian 



given by the following expression 



1^- 



(1) 



+1 



l). 



(3) 



The driven system defined by (01 and Q can be solved 
Tsll and one finds that the ground state does not change 
and the energy flux is zero up to a critical value A = Ac(/i) 
of the driving field. The ground state expectation value 
of the energy flux Je — {Je)/N becomes nonzero for 
|A| > Ac. The resulting phase diagram is depicted in 

Figm 
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FIG. 1: Phase diagram of the driven transverse Ising model 
in the h — \\\ plane where h is the transverse field while A 
is the effective field which drives the fiux of energy. Pairs 
of dual-conjugate points are shown by filled squares, circles, 
and stars; and the line /i = 1 is self-dual as discussed in the 
text. Power- law correlations are present in the nonequilibrium 
phase, Je 7^ 0, and on the Ising critical line in the equilibrium 
phase, Je = (/i = 1, < |A| < 2). 



where Si = and erf {a — x, y, z) denotes the three 
Pauli matrices at sites i = l,2,...,Afofac?=l periodic 
chain (s^_|_j^ = s"), and h is the transverse field in units 
of the Ising coupling (J = 1 is set in the rest of the 
paper) . 

The order parameter of this system is — sf , and 
the ground state of this Hamiltonian changes from being 
disordered (Mr) = for ft, > 1 to ordered {M^) ^ for 
h < 1. The transition point h = he '\s & critical point in 
the universality class of the two-dimensional Ising model. 

In order to drive the above system out of equilibrium, 
one can add a bulk field that drives the energy flux. The 
simplest bulk field is the energy fiux itself [iJl thus, in 
order to produce a nonequilibrium state of Hj one should 
find the ground state of the following Hamiltonian 



h X ^ 



XJp 



(2) 



Here the driving fleld A is again measured in units of J, 
and the current operator J is the sum of local energy flux 



Here we note a hitherto unnoticed property of the 
Hamiltonian Namefy, the duality properties of the 
transverse Ising model 14] have an appropriate general- 
ization to the full nonequilibrium phase diagram of H. 
Indeed, let us denote the action of duality transforma- 



tion by sf 



and define the transformation as it is 



done for the Ising model in a transverse field 



(4) 



It can be verified then that the Hamiltonian H[h, A, {sf}] 
characterized by the two couplings {h, A) transforms 
into an identical Hamiltonian with couplings {h, X)* = 

{h-\-\h) 



H[h,X,{sf}]^hH[h-\-Xh,{sr}]. 



(5) 



Hence the duality transformation leaves the whole |A| = 
Ac curve globally invariant and, furthermore, it leaves 
the h = I line pointwise invariant (examples of dual- 
conjugate points are given in Fig^l. In order to keep 
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the formulas simple, from now on we shall restrict our 
analysis to A > 0, that is to Je > 0. 

The self-dual h = 1 line is expected to display special 
properties. For example, quantities such as Ac/A, or the 
wavevectors q± where the excitation spectrum is gapless, 
are left invariant by the duality transformation. Further- 
more, the functional form of various physical quantities 
(dispersion, energy flux, fluctuations) considerably sim- 
plify on this line and thus the knowledge of selfduality 
helps locating limits where exact calculation can be car- 
ried out. 

Our main goal is to calculate the distribution functions 
P{Mz) and P{Mx) in various regions in the above phase 
diagram. These functions are defined as 



(6) 



where brackets denote the ground-state expectation value 
(note that we have omitted the index a in Pa{■)^ i.e. the 
argument of the function defines which distribution is 
considered) . 

There are two parts to our calculations. The function 
P{Mx) is evaluated numerically by diagonalizing H for 
chains containing up to = 18 — 20 spins, while P{Mz) 
is found analytically. The exact calculation is possible 
because Mz = '}2,q{CqCq — 1/2) is a quadratic form in 
the fermion operators (c^, Cq) in which the Hamiltonian 

H = Hi — XJe is quadratic as well. Thus the calculation 
of the generating function 



G(s) 



(7) 



becomes a problem of evaluating Gaussian integrals. We 
begin with this part of the problem. 



A. Formalism 

The Hamiltonian H can be diagonalized |0| by first in- 
troducing creation-annihilation operators, then employ- 
ing the Jordan- Wigner transformation [I^ to transform 
them into fermion operators (q, c^) and, finally, using 
a Bogoljubov transformation on the Cq and c^^ com- 
ponents of the Fourier transforms of ci-s. The calcu- 
lation of P{Mz) becomes relatively simple if after using 
the Jordan- Wigner transformation one passes to a path- 
integral formulation (see e.g. jl^ for a pedagogical ac- 
count). One then finds that the system is described by 
the following quadratic action 



S[c,c] 




f doj 


fdq 




2tt 



{Cq{uj) C_q{-UJ))B^ 



2 



q^u 



'■-q{-U}) 
Cq{uj) 



(8) 



where the Grassman fields Cg(w) and Cg(a;) are related to 
and Cq correspondingly, while the scattering matrix 



Bq.u, has the following inverse 

{C-q{-Uj)Cq{uj)) (c_g (-w)c_g (-0))) 



{Cq{uj)Cq{uj)) (c^ (w)c_, (-w)) 

Here the correlators are given by 

{Cq{uj)Cq{uj)) = {Cq{uj)Cq{uj))* = 



(9) 



1 

2a: 



I 5 cos g - Ag I 



— AJ 



- I cos g -f Ag 
iu — Aa 



(10) 



and 



{C-q{-Uj)Cq{uj)) = {Cq{uj)C-q{-Uj))* = 



li^ sm q 



2A„ 



iu — Aq iuj — Kq 



(11) 



where Ag and A^ are the dispersion relations for Hi and 
H, respectively 



Kq = -y/l + h'^ + 2hcosq 
Af = ±Ag + —smq. 



(12) 
(13) 



In order to return to the time variables, we must first 
study the two branches A^ of the spectrum. 



B. Spectrum and energy flux 



Flux of energy is pres ent in the system only above a 
critical drive A > Ac [13 (see FigOJ where 



Xcih) 



2 if /i > 1 
I if /i < 1 

h — 



Indeed, it is not hard to see that 

A+(A<Ac)>0 ; A-(A<Ac)<0. 



(14) 



(15) 



Hence the ground state is unchanged with respect to the 
pure Ising model ground state as long as the driving field 
does not exceed Xc{h). This means that all observables 
will assume their Ising model values and no energy cur- 
rent will be flowing through the chain. 

However, for A > Ac one may see that A+ (resp. A~) 
changes sign over the interval I2 = [q-,q+] (resp. I4 = 
[—q+, —q-]). The explicit expression of the wave vectors 
q± is deduced from 



-4± V(A2/i2 -4)(A2 -4) 
cosg± = -7r<<7±<0. 

(16) 

Beyond the critical drive, the excitation spectrum gives 
rise to a ground state which breaks the left-right sym- 
metry and, indeed, it may be verifled |l^ that a nonzero 
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energy flux is present in the chain with the exphcit form 
of the flux is given by 



JE 



47rA2 



(17) 



For later applications we specify the js = J£;(A) function 
in the vicinity of the critical drive, as A ^ A^, 



^^VA^ for/i> 1,A,==2 



3E 



4ir 



3L(A-2) for/i = l,Ac = 2 (ig) 



A-f for/.<l,A, = f 



Besides, as they will naturally arise in the upcoming dis- 
cussion, we further define here intervals Ii = [— tt, (/_], 
^3 — [9+1^9+] ^-iid — ["9-1 ""l which are complemen- 
tary to I2 and I4 in [— tt, tt]. 



and the fields in the exponentials are evaluated at some 
fixed time. In order to evaluate 1)22(1 . we recall the fol- 
lowing result for Grassmann integrals the vacuum 
expectation value of observables of the form e"^? ^I'^i'^i 
can be obtained as: 



n 



/det 



22 



(23) 



As one can see G{s) is a special case of eq. (|23|l and it can 
be evaluated by using the appropriate expressions 1(19(1 or 
for B^i. 

If no current fiows, that is for A < Ac, we find that InG 
(the cumulant generating function) is given by 



In G(s) 



N 



dq 
2tt 



In [(1 - Uq) e" + Uqe 



(24) 



C. Inverse of the scattering matrix in the absence 
or presence of energy flux 

When no current flows through the system, je — 0, 
the equal-time transform of (i3g_i^)~^, denoted by B~^, 
reads: 



I 2 sm q 



2A, 2A, 
■^ + ^cosg+Aq 2^ sing 

2a;; 2A„ 



(19) 



For jE 7^ 0, on the other hand, is given by Eq. p9() 
only for g e /i U /s U 75. If g G J2 U J4, its expression is 
changed to 

h: B-^=( O"), U: B-^=(y^] (20) 
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Having the expressions for B^ ^, we can start the calcu- 
lation of the distribution function P{Mz)- 



III. DISTRIBUTION FUNCTION FOR THE 
TRANSVERSE MAGNETIZATION 

A. Calculation of the generating function and its 
moments 

The generating function ^ of P{Mz) can be expressed 
through fermionic operators as 

G{s) = (e-'**^^) = e^5(e-«i:,4^<.) . (21) 

After normal ordering e"^'^''^' and using the Grassman 
fields, we are left with evaluating the following expression 



G(s) =e^4(e-"^^<<'^«'='), 5 = 1 - e"" (22) 



where 



(h 



5g + 2A,)/(4A,). 



(25) 



One can verify that the normalization condition G(0) = 1 
is satisfied, and one can also recover the well-known result 
found by Pfeuty '17*1 for the magnetization 



5 In G(s) 



ds 



- (M,) = N 



^ dq 



1 



(26) 



As we shall see below, P{Mz) is a Gaussian thus, in 
addition to (Afz), the variance of the transverse magne- 
tization Na^{h) = (M^) — (Mz)^ will characterize the 
distribution. It can be obtained from the second deriva- 
tive of In G(s) as 



"9) 



(27) 



It is interesting to note that the fiuctuations in Mz are 
independent of the magnetic field in the ordered phase 



_( 1/4 



for \h\ < 1 



l/(4/i2) for \h\ > 1 



(28) 



In the presence of nonzero energy flux (A > Ac), the 
generating function is more complicated only because of 
the limits of integration in ((24(1 



lnG,(s) = ^ I gln[(l-n,)e^-t-n,e-^] 



(29) 



/1U/3U75 



Accordingly, the first and second cumulants of the mag- 
netization are given by 



(Af 



z/X 



N 



dq h + cos q 
27r 2AI 



(30) 



/1U/3U/5 



5 



dq 



SIT? q 



I1UI3UI5 



27r 1 + /i2 + 2hcosq 



(31) 



where we have used ea. ^'2b\i to write out the integrands 
exphcitely. Note that we added a A subscript to G, (Mz) 
and in order to indicate that these quantities do de- 
pend on A for A > Ac. 



B. The transverse magnetization distribution is a 
Gaussian 

In order to show that the transverse magnetization dis- 
tribution is a Gaussian, let us consider the n"^ cumulant 
(M")c which is the coefficient in front of (— s)"/n! in the 
expansion of In G{s). The latter coefficient is N times an 
integral of a polynomial in n^, both in the current car- 
rying and current free phases and, furthermore, Uq is a 
nonsingular, strictly positive function of q (note that Ug 
is finite even at h ~ 1). Hence each cumulant depends 
linearly on N . In particular, as we have seen, the vari- 
ance of Mz denoted by iVcr^ and cr^ (in Eq. or (|?7|) ') 
is finite. 

It follows from the linear N dependence of the cumu- 
lant s that 



N 



]yn/2 



7V(2-")/2 



(32) 



thus the above ratio goes to zero for all n > 2 as iV ^ 00. 
We may therefore conclude that the limiting form of the 
distribution function of the transverse magnetization is 
a Gaussian of variance Na'^: 



P{Mz) 



1 



(33) 



The above result applies over the whole phase diagram 
and it is useful to check the numerical procedure em- 
ployed in Section Hvl bv evaluating P{Mz) for finite-size 
systems. As can be seen in Fig|21 there is a convergence 
to the limiting form with increasing N and, furthermore, 
the nearly Gaussian fluctuations of are observed al- 
ready at small sizes {N = 16 — 20). It is remarkable that 
the deviations from the Gaussian are very close to those 
of an N step random walk with a drift determined from 
a correspondence between the left (right) moves and the 
up (down) spins generating the average (M^). 
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FIG. 2: Distribution function P{Mz) for the longitudinal 
magnetization Mz on the critical line h = he — 1, X < 2. 
Results for periodic boundary conditions are displayed. The 
solid line is the asymptotic Gaussian while the dashed line 
is the displacement distribution of a 20-step random walk of 
steplength 1/2 having a drift generating an average displace- 
ment equal to (Mz). 



The above quantities exhibit singular behavior as one 
enters the current-carrying phase. For the magnetization 
we find 



9+ 



5Mz ^ -N 



dq h + cos q 
27r K, 



^(l-/i-2)(A 
dMz ^{ #(A- Ac) 

^/i5/2(A-Ac)3/2 

and for the variance 

h>l,Xc = 2, 5(tI 
5al 



Ac) 



1/2 



for h> 1, Ac 
for /i 1, Ac 
for h < 1, Ac 



Using lfTH|) we find, as ^ 0, 
/i = 1, 5al ~ 



-]E 
.1/2 



(35) 



2 

2 

2 

h 

(36) 



(37) 



(38) 



As can be seen, the variance of the transverse magneti- 
zation is smaller in the current carrying phase than in 
the current free phase. This supports the view that im- 
posing a current stiffens the system, and thus decreases 
fluctuations. 



Recalling the expression of the average and the vari- 
ance (|26I27II we calculate them in the limit of vanishing 
flux (A ^ A+). Let us deflne 

6Ah = {Mz){h, A) - {Mz){h), 5a^ = a^{h, A) - a^{h) 

(34) 



IV. DISTRIBUTION FUNCTION FOR THE 
LONGITUDINAL MAGNETIZATION 

The exact evaluation of P{Mx) appears to be a non- 
trivial task and we have been able to calculate it only 
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numerically for finite size chains. Since the expression 
© for P{Mj.) is a ground-state expectation value, we 
had to find the ground-state wave function and, due to 
the sparsencss of the Hamiltonian matrix, the Lanczos al- 
gorithm 2^] could be used effectively. Since the ground 
state wave function is needed with precision, we were able 
to accomplish this task for chain lengths of up to = 20 
with the results displayed on Figs. I3I5I 



A. Equilibrium distribution 

In the equilibrium system [je = 0), the correlation 
length is infinite only at the critical point. Thus one 
expects P{Mx) to be a Gaussian for ft, > 1, a sum of 
two Gaussians for ft < 1, and a nontrivial distribution 
emerges only at criticality (ft = ftc = !)• This is indeed 
what we observe, apart from the finite size effects showing 
up in nongaussian corrections close to ft = 1. At the 
critical point itself, P{Mx) shows fast convergence to the 
asymptotic form as can be seen on FigI31 where the N > 
16 points appear to have settled on the asymptotic curve. 
This means that the N dependence of P{Mx) is almost all 
in the scaling variable M^/a/ (M^), a remarkable feature 
that has been observed in a series o f eq uilibrium- and 
nonequilibrium critical states HHHHl IH IH . Note 
also that the finite-size effects show up mainly in the large 
Mxl \/ (M^) region. This is in accord with the general 
observation that the large-argument region of the scaling 
function is related to the long-wavelength properties of 
the system. 
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FIG. 3: Scaling function for the distribution P{Mx) of the lon- 
gitudinal magnetization Mx on the critical line h — he — 1, 
A < 2, for periodic boundary conditions. In order to demon- 
strate the smallness of finite-size effects, the small systems 
(A'^ = 8, 10, 12) are displayed by full symbols while the larger 
systems (A^ — 16 — 20) are all shown by a single empty symbol. 

Fig^ shows the critical point scaling functions for var- 
ious boundary conditions (periodic s'^_^_i = s", antiperi- 
odic s'^j^i = — s", and free). One can observe here not 



only the strongly nongaussian character of the distribu- 
tions but also the fact that scaling functions do vary with 
changing the boundary conditions. The boundary condi- 
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FIG. 4: Distribution function P{Mx) for the longitudinal 
magnetization Mx at the critical line h = he — 1, X < 2. 
Results for periodic, free, and antiperiodic boundary condi- 
tions are displayed for system sizes N = 16 — 20. 

tion dep endence of the critical scaling functions is known 
mill El 113 ■ It is also known that the scaling func- 
tions depend on the shape of the system as well. In case 
of the d = 2 Ising model, this means that the scaling 
function depends on the aspect ratio, a, of a rectangular 
sample. Since the transverse Ising model has its origin 
in the transfer matrix of the d = 2 Ising model in an 
anisotropic limit p^ . we believe that the distribution 
functions displayed in Fig0] are equal to the d = 2 criti- 
cal order parameter distributions in the a —^ limit with 
the boundary conditions in the "short" direction being 
the same in the d — 1 and d = 2 systems. Implicit in this 
belief is the assumption that the boundary conditions in 
the "long" direction do not affect the scaling function 
provided a — > 0. 



B. Nonequilibrium distribution 

In the nonequilibrium case {je ^ 0), we find that sim- 
ilarly to Mz, the fluctuations of the longitudinal magne- 
tization are also Gaussian (Fig[SJ|. Remarkably, the 
finite size, non-gaussian corrections are small even near 
the ft = 1, A = 2 point. 

The Gaussian result is somewhat surprising since the 
fiux generates power-law correlations in (s§s^) decaying 
with distance as 1/ ^/n [l^ . Thus one may presume that 
the current carrying states are effectively critical. This 
is not so, however, since the correlations are oscillating 
and the fluctuations (Mj) (given by the integral of the 
correlations) are flnite. Actually (Mj) is decreasing with 
increasing jE (see below) and thus we see here another 
example of fluxes generating power-law correlations but 
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FIG. 5: Distribution function P{Mx) for the longitudinal 
magnetization away from the critical line h — he — 1, 
A < 2. 



Once the correlations are known, the fluctuations can 
be calculated from 



P + OC 

{Ml)=N[l + 2y^C.,{r)]^N drC,(r) (41) 

r>l 



where changing the sum into integral is again allowed 
because of the diverging characteristic length scale . 
Using now the scaling form (|40|l we find that 



{Ml) cx Nk~^'^ oc Nj^""^ . 



(42) 



The above expression demonstrates the decrease of fluc- 
tuations with increasing flux and it also tells us how the 
Gaussian distribution crosses over to the nontrivial shape 
observed at the critical point. 



FINAL REMARKS 



nevertheless making the system more rigid. 

The decrease of (Mj) with increasing can be seen 
in the numerical studies of the finite spin chains. On 
the selfdual line {h = 1), however, this can be shown 
analytically in the limit of vanishing flux [\ —t \c = 2) 
where one finds 

(A^2)cxiVJ//^ (39) 

The derivation of the above result is possible because the 
(/i = 1, A = Ac) point is a critical point with infinite cor- 
relation length. Approaching this point along the /i = 1, 
A — *■ A+ line, one can observe from numerical studies [T^ 
that the wavelength of the oscillation of the corre- 
lation function Cx{r) — (sfsf+^) (x cosKr diverges as 
oc (A — Ac)^^/^. This diverging wavelength allows 
one to take a continuum limit and to establish [s^ to all 
orders in perturbation theory, that the order-parameter 
correlations possess the following scaling form 

C.{r)^^^H^r). (40) 

with the small argument limit of the scaling function ex- 
plicitely given by T{x) = 1 - /2 + x^/l& + 0{x^). 

The derivation of the above results follows the ideas 
[Sll Is^ I33I Is^ used for the calculation of the order- 
parameter fluctuations in the equilibrium transverse Ising 
model. Namely, the correlations are expressed as a Pfaf- 
fian of a block Toeplitz matrix constructed of 2 x 2 
matrices. Then, in the equilibrium case, the analy- 
sis of the Toeplitz matrices in the asymptotic limit of 
r — !■ 00, and h ^ \ with r{h — 1) kept fixed yields 
^x'^i^) r~^/^F{r{h — 1)). A similar asymptotic analy- 
sis in the current-carrying phase, using the scaling limit 
r — > 00, and k ^ (A — Ac)^^^ — > with nr kept finite, 
results in eg . H4()|l . The derivation is rather technical and 
we shall present it, together with the analysis of other 
scaling limits, in a separate publication js^. 



Returning to the problems discussed in the Introduc- 
tion, we can see that the connection between NESS and 
critical states in terms of (universal) distribution func- 
tions is not straighforward. NESS is generated by fluxes, 
and fluxes may or may not generate long-range corre- 
lations. There are numerous examples Q, |^ where 
the fluxes are spatially localized and long-range correla- 
tions do not develop (unless the system is at a special 
point in the parameter space). Clearly, in such cases, 
one cannot hope for a general description to emerge. If 
the fluxes are global, as is the case for the model treated 
in the presentpaper, long-range correlations do emerge 
frequently PjSQ- Even in this case, however, it is far 
from trivial whether these correlations drive the system 
to an effectively critical state or whether they make the 
system more rigid. 

The driven transverse Ising model treated above is an 
example where a global flux of energy generates long- 
range correlations but the resulting state becomes more 
rigid in the sense that the fluctuations are decreased due 
to the presence of the flux. Driven diffusive systems 
1] provide other examples j36j | where the fluctuations 
decrease while the fluxes induce power-law correlations. 
Thus we should conclude that, in general, the power- 
law correlations generated by global fluxes cannot be the 
source of possible universality of nonequilibrium distribu- 
tion functions. It remains, however, an intriguing ques- 
tion whether the weak long-range interactions supported 
by global fluxes can underlie a kind of " weak" universal- 
ity classification of distributions in NESS. 
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